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1. Introduction 

Lattice calculations of light-quark quantities such as quark masses, the kaon bag parameter, 
and semi-leptonic form-factors play a crucial role in extracting the fundamental parameters of the 
Standard Model and in constraining new physics beyond it. More complicated quantities such as 
non-leptonic kaon decays are very sensitive to new physics, but are more difficult to calculate using 
lattice QCD. Nonetheless, great progress has been made on this front in just the last year. 

Almost all of the more straightforward quantities mentioned above have been calculated by 
a number of different groups, and there is excellent agreement between them when we restrict 
ourselves to calculations for which a complete systematic error budget has been provided. Further- 
more, for almost all of the quantities considered here there is more than one result with precision 
comparable to the highest quoted precision for that quantity, so that this agreement is a nontriv- 
ial and important test of lattice methods. This is true of all the light quark masses, including the 
strange quark mass, where it appears that the discrepancies of the past have now been largely 
resolved. (See, e.g. Ref. [|l|] for a compilation of older quark mass results.) 

In this talk I cover light quark masses (u, d, and s), the K — > nlv form factor and the de- 
termination of \V US \, the kaon bag parameter B K , kaon bag parameters relevant for beyond the 
Standard Model physics, and K ^ Tin matrix elements. There are fully documented results from 
more than one group for most of these; when this is the case world averages are presented. The 
averages presented here is work done in collaboration with Enrico Lunghi and Ruth Van de Wa- 
ter, where we update and extend to new quantities the averages presented in reference The 
latest averages (including quantities outside the scope of this review) can be found on our website 
www.latticeaverages.org. After reviewing the procedure we adopted for the averages presented in 
this talk, I present status updates of the various quantities in turn, including the new averages. 

Related reviews presented at this conference include a review of the light pseudoscalar decay 
constants by Christian Holbling [^J, and a review of the impact on phenomenology of many of the 
above quantities by Chris Sachrajda [Q|. 

2. Lattice averages 

Our main criteria for including a result in the world averages quoted here is that the calculation 
include a complete error budget, including all sources of systematic error, and also that the work be 
documented in either a publication or proceedings. Thus, we do not include numbers that appear 
only in the slides of a conference talk. For the most part, we restrict ourselves to including only 
Nf = 2 + 1 flavor lattice results in the averages because Nf = and 2 flavor calculations typically 
do not quote an error due to quenching, as this is notoriously difficult to estimate reliably. Since we 
only include results where estimates of all relevant systematic errors have been made, we exclude 
most of the Nf = 2 flavor results from the averages, though I show the most recent state-of-the-art 
two flavor calculations in the comparison plots. We make an exception in the averages for the 
Nf = 2 calculation of the K — > nlv form factor at zero-recoil, /+(0), since chiral perturbation 
theory power counting can give a reasonable estimate for the size of the effect of quenching the 
strange quark for this quantity, and this error was included as part of a complete systematic error 
budget by the ETM Collaboration in ref 
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There is an additional small systematic error in all of the Nf = 2+1 flavor lattice results due to 
the neglect of the dynamical charm quark from the simulations. This error is typically not included 
in current error budgets, though the effect is believed to be small. Because the charm quark is 
relatively heavy compared to Aqcd, the size of the effects due to neglecting the charm quark in 
the sea can be estimated using heavy-quark effective theory (HQET). The leading corrections are 
of order a s [A/(2m c )], where the two is a combinatoric factor appearing generically in HQET. 
Using reasonable estimates for these quantities (e.g., a s = 0.33, A = 500 MeV, m c = 1.2 GeV) 
suggests that the size of charm loop effects is around 1%, although this effect is likely to be further 
suppressed in SU (3)-breaking ratios such as JkI f%- Simulations with Nf = 2 + 1 + 1 flavors are 
underway by the ETM Collaboration [||] and the MILC Collaboration [f7|] to address this issue. 
Note that this HQET power counting estimate cannot be applied to the strange quark, and that the 
effects due to quenching the strange quark are more difficult to quantify without a direct comparison 
between Nf = 2 and 2+1 flavor simulations. Such a direct comparison is possible, since the ETM 
Collaboration has now calculated many quantities using Nf = 2 flavors with all systematic errors 
under control except that due to quenching the strange quark. Good agreement is found with 2+1 
flavor calculations at the quoted level of precision, suggesting that the effect of quenching the 
strange quark is small. Nonetheless, the precision on the 2 + 1 flavor averages is typically better 
than that of the 2 flavor calculations and is of the size where one might expect dynamical strange 
quark effects to become visible. Thus, to be consistent with our criteria of including only results 
with estimates of all the relevant systematic errors, we do not include Nf = 2 flavor calculations 
in our averages, expect for the special case of /+(0), as mentioned above and discussed further in 
Section V. 

In order to avoid underestimating the errors in the averages, we need to take into account 
correlations between the different calculations. Although the full correlation matrices between 
the different calculations do not exist and would be difficult to construct, we can still account for 
correlations using the following conservative assumptions: when calculations use the same gauge 
field ensembles for the same quantity, we assume the statistical errors are 100% correlated. When 
a systematic error is at all correlated between different calculations, we assume the correlation in 
this error between calculations is 100%. This treatment of systematic correlations is conservative 
in that it will lead to somewhat of an overestimate of the total error in the average, but we feel it is 
a reasonable assumption given the information available. Finally, we also adopt the Particle Data 
Group (PDG) prescription for combining several measurements whose spread is wider than what 
one would expect given the quoted errors. The error on such an average is rescaled by the square 
root of the minimum of the % 2 / per degree of freedom 

A snapshot of the simulation parameters is presented in Table [I] to give a sense of the range of 
lattice actions, lattice spacings, and light pion masses being used in the most recent results. For the 
staggered simulations, I have quoted the root-mean-square pion mass for the minimum pion mass 
in the sea, while for the valence sector I have quoted the taste-goldstone pion mass as the minimum 
pion mass. This gives a sense of how light the valence and sea pion masses are compared to other 
simulations and is based on a more detailed look at staggered chiral perturbation theory formulas 
fll7| , [l8| , |59| ] for the quantities considered in this review that were computed using staggered quarks. 
For simulations performed with domain wall valence quarks on a staggered sea, I again quote the 
root-mean-square pion mass for the minimum pion mass in the sea, while for the valence sector 
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Table 1: Snapshot of parameter values being used in numerical simulations reviewed in this talk. The first 
column is the group name; the second is the number of dynamical fermions; the third is the action; the 
fourth is the approximate range of lattice spacings; the fifth is the dimensionless product m K L, which gives a 
measure of expected finite-size effects; and the last column is the approximate minimal pion mass simulated, 
in both the sea and valence sectors. 



Group 


N f 


action 


a(fm) 


m n L m™ m (MeV) 










sea/val 



ETMC 



MILC [10] 



RBC/UKQCD [11 



JLQCD [12] 



PACS-CS [13] 



BMW [14] 



ALV[70] 



HPQCD [16] 



2 Twisted Mass 

2+1 (Asqtad) staggered 

2+1 Domain Wall 

2+1 Overlap 

2+1 Clover 

2+1 Clover 

2+1 DW on MILC 

2+1 HISQ on MILC 



0.05-0.10 fm > 1 280/280 

0.045-0.12 fm >4 250/180 

0.085-0.11 fm >4 290/210 

0.11 fm >2.7 310/310 

0.09 fm > 2.0 140/140 

0.065-0.125 fm >4 120/120 

0.06-0.12 fm > 3.5 250/210 

0.045-0.15 fm > 3.7 360/310 



I quote the mass of the valence pion made of two domain wall quarks. This table reflects the 
parameters used in ongoing simulations, not necessarily the parameters appearing in the averages 
for all quantities presented below. 

3. Quark masses 

The quark masses are fundamental parameters of the Standard Model, but due to confinement 
do not appear as free particles in nature. One must tune the lattice quark masses so that one 
reproduces the experimentally known hadron spectrum, with three experimental inputs needed to 
determine three quark masses (and a fourth experimental input to fix the scale). The bare lattice 
quark masses are then known in the regularization scheme defined by the lattice action of any 
particular calculation, but to compare to results using other lattice actions or to be used as input 
in continuum calculations, must be converted to a standard continuum regularization scheme like 
MS. The matching factors needed to convert to the MS scheme are short distance parameters, and 
can be computed perturbatively. Direct perturbative matching between a lattice scheme and the MS 
scheme is a technical challenge, especially at two-loop order or beyond. The convergence of lattice 
perturbation theory is typically poor unless tadpole-improvement and a renormalized coupling are 



used [20]. The convergence also tends to be better behaved when improved actions with fat gauge 



links are used [2lp, but this increases the difficulty of going to higher order, and few results beyond 
one-loop are known. Non-pertubative matching to a regularization independent continuum scheme 
can be performed [p2|], and then the matching to MS can be done in continuum perturbation theory, 
where it is easier to go to higher orders. It is also possible to renormalize the masses completely 
non-perturbatively using the Schrodinger functional scheme [23, Qj. 



The latest results for the strange quark mass and for the averaged up and down quark mass 
m uc i are shown in Figure |T} The notation in these and all subsequent average plots is as follows: for 
quantities that have been presented with full error budgets in a paper or proceedings, the data point 



4 



Light quark physics 



Jack Laiho 




MILC '09 
HPQCD '10 
RBC/KEK/Nagoya '10 
RBC/UKQCD '10 
BMW '10 
ALV '09 
MILC '10 
PACS-CS '10 
ETMC'10(2 flavor) 



80 90 100 

MS(2GdV) 

m (MeV) 



110 



120 



I O- 

I © I 



I — — I 



O MILC '09 

• HPQCD '10 

• RBC/KEK/Nagoya ' 10 
9 RBC/UKQCD '10 

• BMW 10 
O ALV '09 

© PACS-CS '10 

• MILC '10 

9 ETMC '10 (2 flavor) 



2.5 



35 4 
MS(2 GeV) 



4.5 



5.5 



(MeV) 



Figure 1: Left panel: The strange quark mass m^ 5 (2GeV) from many groups, including the world average. 
See the text for further details on color coding, etc. Right panel: The results for the average of the u and 
d quark masses »r™?(2GeV) from many groups with world average. Results in both figures are taken from 

Ref s. is ra m pi pi |§ m m m Hi . 



is shown in green, and is included in the average. Quantities that are not included in the average 
are shown in red. This can be because the quantity has not yet been fully documented, or because 
potentially significant systematic errors have not yet been estimated. If the situation is the later, the 
quantity is distinguished by a dashed line in the error bar. A slight exception to this rule is that two 
flavor results are displayed without a dashed line, to emphasize that all errors except that due to 
quenching the strange quark have been included in the error estimate, but they are not included in 
the averages in most cases as discussed above and are labeled explicitly. 

The agreement between the different results for quark masses is striking, especially given 
the long history of disagreement between various m s determinations. The only noticeable outlier 
in Figure |] is the PACS-CS result [31], though this result does not yet include a full systematic 
error budget by the admission of the authors, so that the significance of the difference cannot 
yet be assessed. The PACS-CS result uses a fully non-perturbative renormalization of Z m from 
the Schrodinger functional method, and values of the quark masses down to the physical masses. 
These simulations still involve rather small volumes {m n L ~ 2) and there is so far only a single 
lattice spacing. The quoted systematic errors include statistical errors, a systematic error due to 
reweighting to the physical quark masses, and a complete systematic error for Z m . The error budget 
does not yet include errors due to discretization effects and finite volume effects, but calculations 
to address the remaining systematic errors are in progress 

A few results in the quark mass determinations stand out due to their especially small errors. 
The quark mass determinations of HPQCD use a novel method and have the smallest quoted er- 
rors [|6|j. The HPQCD Collaboration make use of their determination of the charm quark mass 



using current-current correlators [ }Z6| p4| ], and the ratio of charm to strange quark masses [16] us- 
ing a highly improved staggered quark (HISQ) action [531] to simulate a relativistic charm quark 
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0.000 0.005 0.010 0.015 0.020 

a 2 (in fm 2 ) 

Figure 2: Continuum extrapolation of ratio m c /m s from HPQCD ll 



so that the charm and strange quark can be treated in the same formulation. The ratio m c /m s , 
which receives no renormalization if both quarks are simulated in the same relativistic formalism, 
allows HPQCD to propagate their precise charm quark mass determination [ |26| ] (using the HISQ 
formalism and four-loop continuum perturbation theory [[36]]) down to the light quark masses. The 
continuum extrapolation of this quantity from HPQCD is shown in Fig. ^[ The m c j m s ratio has also 
been calculated by the ETM Collaboration using two dynamical flavors [^2|], and it is in good agree- 
ment with the result of HPQCD. The BMW Collaboration uses smeared clover quarks down to the 
physical light quark masses and the Rome-Southampton non-perturbative renormalization (NPR) 
method [|22|] and several lattice spacings down to ~ 0.05 fm to control the perturbative match- 



ing [28]. The RBC/UKQCD Collaborations use domain wall quarks and the Rome-Southampton 
NPR method, but they match to a non-exceptional momentum scheme [11, 37], reducing the con- 
tamination from infrared effects and significantly reducing the size of the perturbative corrections 
compared to schemes that involve exceptional momenta. 

Although the average of the u and d quark masses can be obtained in a fairly straightforward 
way from lattice simulations, the difference between the u and d quark masses requires an under- 
standing of electromagnetic effects. These have been incorporated into lattice calculations using 
different approaches. The MILC results for m u and use the differences between the masses of 
charged and neutral pions and kaons, along with continuum estimates of the violation of Dashen's 
theorem to estimate the difference between the u and d quark masses [38]. HPQCD makes use of 
the MILC values of the ratios between the u and d quark masses, thus using the same input for 
electromagnetic contributions as MILC. RBC/KEK/Nagoya include (quenched) QED explicitly in 
the lattice calculations, thus bypassing the need for continuum information on electromagnetic ef- 
fects [27]. They make use of two volumes to study finite volume effects, and they use SU (2) heavy 
kaon chiral perturbation theory [ ^9| ] to perform the chiral extrapolation, including electromagnetic 
effects. BMW obtain information on the electromagnetic corrections from dispersive studies of 
v\ — y 371 decays [14]. Again, the agreement is impressive, as can be seen in Figure [| and the errors 
on the world averages presented here for the u and d quark masses are at the 5% and 2% level, 
respectively. Note that the u quark mass is over 20a from zero, so that the vanishing of this mass 
appears to be ruled out as a solution to the strong CP problem [[f(], 41, 43, 42]. 
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Figure 3: Left panel: The u quark mass m„ (2GeV) from different groups, including the world average. 
Right panel: The d quark mass m^f s (2GeV) from different groups with world average. This figure quotes 
results from Refs. [12. 



Figure |] shows the quark mass ratios m s /m uc i and As noted above, these ratios are 

interesting to examine because the renormalization factor needed to quote the quark masses in a 
common scheme cancels. For m s /m Ui i there is, once again, impressive agreement between the var- 



ious determinations, in particular those of MILC [g5p, HPQCD [Jig], and BMW where the 
total errors from each group are at the sub-percent level. This precision agreement is noteworthy in 
part because it is not possible to get much guidance from experiment on the value of this ratio, un- 
like say the light hadron spectrum or pion and kaon decay constants, and no other non-perturbative 
methods quote errors that are competitive with the precision of the lattice results. Leading order 



chiral perturbation theory (#PT) predicts a value around 26.0 [44], but it is difficult to improve this 
result systematically to the same precision using #PT because of uncertainties in the low energy 
constants appearing at higher order in the chiral expansion. 

4. K — > 7t£v semileptonic decay 

The semileptonic K — > nlv decay can be used to obtain the CKM matrix element V us from the 



experimental branching fraction using [ 45 1 
„ Glm 



KB 



T^C|5 w (K|/ + ) 2 /k(1 + 5«+5 5 ^ (2) ) 2 , (4.1) 



where Sew = 1-0232(3) is the short-distance electroweak correction, Ck is a Clebsch-Gordan co- 
efficeint, /+ (0) is the form factor at zero momentum transfer, and Iki is a phase-space integral that 
is sensitive to the momentum dependence of the form factors. The quantities and are 
long-distance EM corrections and isospin corrections, respectively. The value 

|Vj/ + (0)= 0.2163(5) (4.2) 
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Figure 4: Left panel: The ratio of quark masses m s /m m j from different groups, including the world average. 
Right panel: The ratio of quark masses m u /m,i from different groups with world average. Results in these 
figures are taken from Refs. [|||l6[|6[|7}|llJ|8||9}|o||Tj|32|. 
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Figure 5: Results for /+(0) including our average. The results are quoted from Refs. [|] fl9| ], 



has been determined from experimental measurements of K — > %i\ decays and non-lattice theory 



for the other inputs to Eq. ( ft. 1[ ) [}45[]. The non-perturbative information is encoded in the form 
factor /+(0), and once this is known from lattice QCD, a value for V us can be determined. The 
value of /+(0) is already rather well constrained by SU (3) chiral perturbation theory, an expansion 
in powers of m 2 K /{%K 2 f 2 ) [[?6|]. One can write /+(0) = 1 +/2 +Ja + where the first term is equal 
to one due to current conservation in the SU (3) limit. The correction fj, does not contain any new 



unknown low energy constants, as required by the Ademollo-Gatto theorem [47], and is predicted 
in terms of pion and kaon masses and the pion decay constant to be f% = —0.0226. We need to 
know /4, if we are to do better, but this requires the determination of new higher order unknown 



low energy constants. The value for f\ was estimated by Leutwyler and Roos in 1984 [48] using a 
quark model; they obtained f 4 = -0.016(8), which gives /+(0) = 0.961(8). 

There are only two lattice calculations of /+(0) that have complete systematic error budgets, 
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one from ETMC from 2009 |Q] and a more recent update in 2010 [49] of an older calculation from 
RBC/UKQCD [5C], with the average shown in Figure |5[ As mentioned in Section^, we include the 
ETMC result in our world average even though it is with Nf = 2 flavors because there is a credible 
estimate of the systematic error due to quenching the strange quark. The ETMC result has already 
been reviewed by Vittorio Lubicz [ p3| ] in his lattice review talk from last year, but it is worth 
recalling the treatment of the error due to quenching the strange quark for /+(0). In the chiral 
effective theory the Nf dependence can be incorporated; this modifies the chiral logarithms and 
changes the low energy constants from their Nf = 3 to Nf ■ = 2 values. Although it is straightforward 
to use the known chiral logarithm expressions to adjust for unphysical quenching effects, the Nf 
dependence of the low energy constants is not known. For / + (0) the effect of using Nf = 2 sea 
quarks can be corrected for in the fa term exactly, since there are no new low energy constants 
appearing to this order, and the strange sea quark dependence is a known expression involving 
chiral logarithms. The fa term requires the knowledge of new low energy constants, and these 
will take their Nf = 2 values in the ETMC calculation. The remaining error in /+ (0) due to the 
difference between fa in the Nf = 2 and 3 flavor theories is estimated by looking at the difference 
between fa in the Nf = 2 and Nf = calculations. The full difference between fa in the 2 flavor 
theory and the flavor (completely quenched) theory is taken as the error due to quenching the 
strange quark for this quantity [|J. 

The RBC/UKQCD result for / + (0) fg|] is an update of their previous result [f>0p from 2007. 
This new calculation improves upon the previous one by using twisted boundary conditions 52] 



to remove the need to interpolate in q 2 to get the form factor at zero-recoil. They also apply a 
different choice of chiral extrapolation, where analytic next-to-next-to leading order (NNLO) terms 
that do not obey the Ademollo-Gatto theorem are included. Such terms are possible if f n is used 
in the NLO expression instead of fa, the decay constant in the SU (3) limit. This is because when 
reordering the series to use fa in the NLO expression, analytic terms that do not respect the mass 
interchange symmetry can appear at NNLO. RBC/UKQCD quote / + (0) = 0.9599 (34) (+^)( 14), 
where the first error is statistical, the second is due to the chiral extrapolation, and the third is an 
estimate of discretization errors [p9[|. 

There were progress updates for / + (0) from the FNAL/MILC Collaborations using staggered 



quarks Q54{ ] and from the JLQCD Collaboration using overlap quarks [p5Q - The FNAL/MILC 



calculation is using the method developed by HPQCD for D — > K£v [ pq ] to get a result for / + (0) 
using 



/+(0) =/o(0) = -{ !<rc|S|*V =0 (4.3) 

m K~ m n 

so that no renormalization is required. This avoids the use of non-local vector currents and does not 
require multiple three-point correlators to form various double ratios. Twisted boundary conditions 
are used to calculate fa at the q = point. The disadvantage of this method is that one cannot 
obtain fa(q 2 ) for q 2 / 0, but this is still sufficient to determine \V US \, since the shape dependence of 
the form factor is usually taken from experiment. 

The JLQCD calculation was done at a single lattice spacing and a somewhat small (1.7 fm) 
volume. The q 2 dependence was modeled using various functional forms to interpolate to q = 
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Figure 6: Left panel: Extrapolation in light quark mass for B% from RBC/UKQCD ||5 
Comparison of Bk in different renormalization schemes from RBC/UKQCD 



Right panel: 



0, and the shape dependence was in reasonable agreement with experiment. Work with twisted 



boundary conditions and larger volumes is in progress [p5 1. 
5. Kaon mixing 

The constraint on the unitarity triangle coming from kaon mixing can be expressed as 

\e K \ =C e K e B K A 2 lf{-r ll S (x c )(\ - X 2 /2) + T] 3 S Q (x c ,x t ) + r/ 2 5 () (x ; )A 2 A 2 (l -p)}, (5.1) 
where C e is a collection of experimentally determined parameters, K e represents long-distance 



contributions and a correction due to the fact that the quantity <j) e / 45 degrees ||57|], Bk is the kaon 
bag parameter, the ty-So are perturbative coefficients, and A, A, p, tj are CKM matrix elements in 
Wolfenstein parameterization. The experimental determination of |ejg-| leads to a constraint on the 
unitarity triangle in the shape of a hyperbolic band in the p-77 plane. The main non-perturbative 
input needed from the lattice to implement this constraint on the CKM unitarity triangle is the kaon 
bag parameter Bk- 

There have been three recent updates on this quantity: from the RBC/UKQCD Collaborations 



[(580, from the SBW Collaboration [J59p, and from ETMC [jSOp. The RBC/UKQCD Collaborations 
have updated their eariler result [^TJ using domain-wall quarks with a number of improvements, 
including a second lattice spacing to allow a continuum limit to be taken, and the use of various non- 
perturbative renormalization schemes with non-exceptional momentum to perform the matching to 



the continuum [p2[]. The RBC/UKQCD results have also modified their approach to the chiral 
extrapolation. For their central value, they average the result using SU (2) heavy-kaon #PT to 
perform the extrapolation in light quark mass and the result using a simple linear extrapolation. 
This approach was motivated by the absence of detectable curvature in their data, and the tendency 
of their SU (2) fit to undershoot f % . Figure || shows the light-quark mass extrapolation for Bk for 
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both the linear and the St/ (2) fits. RBC/UKQCD used multiple RI-SMOM schemes with non- 
exceptional momenta to determine the matching factor for Bk ■ Figure ^ shows the comparison of 
results between the various RI-SMOM schemes, as well as results from the RI/MOM scheme. The 
systematic error associated with the NPR matching is also shown in Fig ^. The final result quoted 
for Bk is 0.749(7)(21)(3)(15), where the errors are statistical, chiral extrapolation, finite volume, 



and renormalization [58]. 



The SBW Collaboration has adopted a mixed-action approach, using HYP-smeared staggered 
quarks [ ]63| ] on the MILC asqtad ensembles, with four lattice spacings down to 0.045 fm [59, 



The matching is done to one-loop order in lattice perturbation theory J65|]. The chiral and contin- 
uum extrapolation is done using SU (2) staggered chiral perturbation theory [|6|. The SU (2) for- 
mulation has much simpler expressions than the SU (3) case, where many new parameters specific 



to the staggered formalism enter Q18|]. SBW quote as their main result a value Bk = 0.724(12) (43), 
where the errors are statistical and the sum of systematic errors in quadrature [|9j. The domi- 
nant error is currently due to the one-loop perturbative matching. Non-perturbative matching is in 
progress. 

The ETM Collaboration has a result for Bk using a mixed action [|7|] with three lattice spac- 



ings down to 0.07 fm [pO[]. The valence action is Osterwalder-Seiler [|58j] and the sea sector is that 
of the ETMC N f = 2 twisted mass ensembles [g]. They quote a value of B K = 0.733(29)(16), 
where the first error contains statistical, chiral extrapolation/fit, and matching errors, while the sec- 
ond contains an error due to different assumptions of &(a 2 p 2 ) dependence in the RI-MOM scheme 



matching factor [60]. A calculation by ETMC using 2+1+1 flavors is in progress. 

The world average for Bk is shown in Fig. [7|. All of the results are in good agreement, which 
is impressive given the different discretizations and methods employed in the various calculations. 
The effect of quenching the strange quark appears to be rather small, as we can see by comparing 
the ETMC result with that of the 2+1 flavor average. 

The ETM Collaboration has presented the first preliminary unquenched (two-flavor) results 
for four-quark operators that contribute to kaon mixing in the presence of new physics [[71]]. The 
effective Hamiltonian relevant for beyond the Standard Model physics contains several new four- 
quark operators, in addition to the one associated with Bk- Only four new matrix elements are 
required, however, since for the new operators only the parity-even parts are needed because the 
strong interaction conserves parity. 

The bag parameters associated with the new operators are defined by [ |7~ij ] 

(^\^)\K°)=B K (n)^m 2 K fl (5.2) 



,-ftO 



(T'\0 i (nW)=C i B i (n) 



m K f K 



m s {n) + m d {n) 



(5.3) 



where Q = {—5/3, 1 /3,2,2/3}, i = 2, ...,5, when we take the basis of operators used in Ref. [|7T|]. 
The matrix element of ff\ is just the Standard Model contribution to kaon mixing. Note that the 
other operators do not vanish in the chiral limit. 

The ETMC calculation includes three lattice spacings, fairly light pion masses (down to ~ 280 
MeV) and non-perturbative renormalization using the RI-MOM scheme. They find that the chiral 
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A 

Figure 7: Results for Bk including the world average. Results are taken from Refs. [|9|, [59[ ||, | 



Table 2: Preliminary values of BSM kaon bag parameters in the MS scheme with two dynamical flavors 
from ETMC [^TJ. Errors are statistical only. 



i 


Bi 


2 


0.56(0.04) 


3 


1.43(0.13) 


4 


0.76(0.06) 


5 


0.63(0.09) 



extrapolation is not very sensitive to the choice of fit function, and we show the results for the B, 
parameters from the quadratic fit of Ref. [|7l|] in Table [| The errors are statistical only. 

6. 



Lattice calculations of non-leptonic K^%% decays are challenging because of the two4iadron 
final state, but they are important for phenomenology. Lattice calculations of K — > %% matrix 
elements have the potential to give us a first principles determination of the decades old AI = 
1/2 rule, and would finally allow us to use the experimental measurement of e' /e as a precision 



constraint on the Standard Model [72]. The Standard Model prediction for e'/e is 

"Im(A 2 ) 



Ref^ 



CO 



V2\e K \ 



Im(A ) 

Re(A 2 ) Re(A ) 



(6.1) 



where Aq and A 2 are the amplitudes for K — )• %% decays into definite isospin states, and real and 
imaginary refer only to the part of the amplitude that becomes complex due to the presence of the 
weak phase. The smallness of the parameter CO = Re(A 2 ) /Re(Ao) ~ 0.05 is a manifestation of the 
A/ = 1/2 rule. 
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K — >■ nn matrix elements are difficult to calculate on the lattice because the Maiani-Testa no- 



go theorem [|73j tells us that we cannot extract physical matrix elements from Euclidean correlation 
functions with multi-hadronic final states. Due to the restriction of working in Euclidean time, the 
most straightforward lattice implementation of calculating K — >• Tin matrix elements only works 
if the final state pions are at rest, or at some other set of unphysical kinematics. Two general 
strategies have emerged for getting around this problem. One strategy is to construct K — > nil 
matrix elements indirectly using the low energy constants (LEC's) of chiral perturbation theory as 



determined from simpler lattice matrix elements such as K — > and K — > % [74]. It was shown in 



Refs. [ 75 , 76, 77 , 78]] that all LEC's through next-to-leading order could be obtained from relatively 
simple lattice quantities. However, this method has the disadvantage that the convergence of SU (3) 
chiral perturbation theory at the physical kaon mass is slow, and it is not clear whether K — > Tin 
matrix elements can be computed in this way to a useful precision [f7^, |80| ]. 

A method for calculating K ^ titi matrix elements directly at physical kinematics was in- 
troduced by Lellouch and Liischer [|8T|]. The Lellouch-Liischer method exploits the finite lattice 
volume to obtain the matrix elements directly by tuning the volume so that the first excited state 
of the two pion state matches the kaon mass. The direct method is straightforward to implement, 
though it is computationally demanding because it requires large lattice volumes (~ 6 fm) and phys- 
ical light quark masses. Improvements to the method have been introduced so that the non-zero 



momentum pion state becomes the ground state and smaller volumes can be used Q82 , S3] , |84 , p5p . 
The RBC/UKQCD Collaborations have made significant progress using the direct method, with 
preliminary results at nearly the physical quark masses and physical kinematics for the M = 3/2 
decay channel [ |86[ ] ; this is discussed below. 

A technique for improving upon the indirect method was presented by the author and Van de 
Water, where we exploited the fact that one can simulate K ^ Tin matrix elements with the pions 
produced at rest [p0||. If one takes the pion mass to be l/2mx, this amplitude can be computed 
directly, since it is not forbidden by the Maiani-Testa theorem, a fact known for quite some time 
J87j|. Thus, we tune the light quark masses such that mx = m^ hys and m n = /n^ hys /2. One can then 
correct for the unphysical kinematics using fixed order SU (3) #PT, where the low energy constants 
can be obtained from simpler quantities, like K — > %. Because the kaon is tuned to its physical 
value, the terms involving kaons (etas) are correct (nearly correct) to all orders in the SU (3) chiral 
expansion. Thus the 10-30% precision of NLO SU (3) #PT now appears in a small correction 
factor, rather than the entire amplitude. This can be tested for known quantities like fx and f n , 
results of which are shown in Fig. || These plots illustrate the fact that the NLO corrections to 
f n are below 10% for this method, and the corrections to fx are below 5%. Also, the one-loop 
corrections account for most of the difference between the 2m n = mx values and the accepted 
values. 

The value for Re(A2) at physical kinematics from JL and Van de Water is calculated similarly, 
where the value at the 2m n = mx point is corrected using the leading order #PT expression, 

(^^|0( 27 . 1 )'( 3 / 2 )|^) LO = ^(ml-ml), (6.2) 
where fo and Bq are the pion decay constant and Bx in the SU (3) chiral limits, respectively. The 
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Figure 8: Demonstration of the size of NLO SU(3) corrections to quantities evaluated at 2»% = m^ hys for 
f K (left) and f% (right). Errors on the circular points are statistical only. Uncertainties in the vertical f K and 
/k error bands include both statistical and systematic errors Q8( 
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Figure 9: The correction to Re(A2) at 2m n — nv K y ' using leading order SU (3) %PT compared to experiment 

[pot. 



leading order correction factor is then 



cLO 



[{m K /2f -m\\l\m\-{m K l2f 



(6.3) 



which is only 23%. This correction is shown in Fig. ^[ along with the experimental value. This 
corresponds to a value of Re(A2) = 1.568(86) x 10~ 8 GeV, where the error is statistical only [pO|]. 
A preliminary estimate of the systematic errors in this approach are shown in Table ||| The error is 
dominated by the 12% #PT truncation error, though this may improve when the correction factor 
in Eq. ( |(0| ) is known to one-loop in #PT. 

The calculation of K — > %% decays from RBC/UKQCD uses the direct Lellouch-Luscher ap- 
proach, and they have made significant progress, including a preliminary result for matrix elements 
in the AI = 3/2 channel with close to 10% errors These calculations are done using the new 
Dislocation Suppressing Determinant Ratio (DSDR) domain-wall quark ensembles being gener- 
ated by RBC/UKQCD \M\ with volumes of 32 3 x 64 with L s = 32 at an inverse lattice spacing 
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Table 3: Estimated total error budget for Re (A2) from JL and Van de Water [p0[]. Each source of uncertainty 
is given as a percentage. 



uncertainty Re(A 2 ) 

statistics 4.7% 

#PT truncation error 12% 

uncertainty in leading-order LECs 4% 

discretization errors 4% 

finite volume errors few percent 

renormalization factor 3.4% 

scale and quark-mass uncertainties 3% 

Wilson coefficients few percent 

total less than 20% 



Table 4: Estimated total error budget for Re(A2) from RBC/UKQCD [|8(|. Each source of uncertainty is 
given as a percentage. 



uncertainty Re(A2) 

statistics 5.8% 

scaling violations 8.5% 

finite volume effects 7% 

partial quenching 2% 

pion phase shift 2% 

meson masses and 2-pion energies 1 .2% 

total 11% 



of around 1 .4 GeV. This corresponds to a spatial box size of around 4.5 fm. The lightest unitary 
pion mass is 180 MeV, and a lighter valence pion with mass around 140 MeV is used for the cen- 
tral value. The main errors contributing to the RBC/UKQCD calculation of Re(A2) are given in 
Table ||. The largest error is the estimate of scaling violations due to the use of somewhat coarse 
lattices at a single lattice spacing and the fact that K ^ Tin matrix elements scale as the lattice 
spacing cubed. The RBC/UKQCD result for Im(A2) is expected to have a similar error once the 
nonperturbative renormalization is completed. 

The calculation of the AI = 1/2 rule is more difficult for a number of reasons. One is the 
presence of power divergent contributions arising from mixing with lower dimensional operators. 
This problem has been addressed by the use of chiral fermions, where the operator subtraction 
is straightforward [89, 50 1. Another problem is the presence of enhanced finite- volume effects 
that afflict the calculation when the light valence quark masses are not the same as in the sea 



>1 , 92]. This was an especially serious problem for quenched attempts to calculate AI = 1/2 kaon 
matrix elements [p3||, but is under control when sea quarks of the correct mass are included in the 
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calculation. Another difficulty is the appearance of disconnected quark flow diagrams, leading to 
the need for very high statistics. The contractions at the level of quark flow are shown in Fig. 10 for 
the AI = 1/2 channel. The red circle is the insertion of the four-quark operator. Additional diagrams 
with a quark current insertion are not shown, but are needed to perform the power divergent operator 
subtraction. Figure |l0|(d) shows the disconnected diagram that is problematic due to the need 
for high statistics. This problem will likely be solved by more computing and better inversion 
algorithms. 

Results for K ^ %% correlation functions in the AI = 1/2 channel from RBC/UKQCD were 
presented by Qi Liu [34], and they are shown in Fig. O. These correlators were obtained from 2+1 
flavor domain wall ensembles with Iwasaki gauge action generated by RBC/UKQCD with volume 
16 3 x 32 and L s = 16, with pion masses around 420 MeV. Propagator inversions were performed on 
each time slice for 400 configurations. The operator Q2 gives the dominant contribution to Re(Ao) 
at renormalization scales typical of lattice calculations, while gives the dominant contribution to 
Im(Ao). The blue (open) circles show the contribution to the correlator from the diagrams (a), (b), 
and (c) in Fig. [T(| while the red (filled) circles include all four diagrams, including the disconnected 
diagram Fig.|0|(d). The comparison between open and closed circles in Fig. fil] shows the large 
statistical errors introduced by the disconnected diagram. Table || also illustrates this point with 
entries for Re(Ao) and Im(Ao) with and without the disconnected diagram. The error on Re(Ao) 
is about 25%, while the error on Im(Ao) indicates that more statistics are needed to be sure of 
a signal. It was reported at the conference by Liu that for non-zero momentum the AI = 1/2 
correlators barely had a signal, even without the disconnected diagrams. Again, improvements are 
expected by going to larger lattices, as well as bigger machines and better inversion algorithms. 

Finally, I observe that the indirect method for K — > %% discussed above is not likely to match 
the precision possible with the direct method in the AI = 3/2 channel, but could be useful for the 
AI = 1/2 channel, given the significant amount of computing needed to reach non-zero momentum 
at the physical kinematics in this channel. 
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Figure 11: Left panel: K — > %% matrix elements at zero momentum for £>2 with and without disconnected 
diagrams from RBC/UKQCD [53]. Z?z'g/z? panel: K — > tttt matrix elements at zero momentum for £>6 with 



and without disconnected diagrams from RBC/UKQCD [94-]. 



Re (A o) no discon 


Re(A () ) 


Im(Ao)no discon 


Im(A ) 


38.7(2.1) x 1(T 8 


30(8) x 10~ 8 


-63.1(5.3) x lfr 12 


-29(22) x 10~ 12 



Table 5: Results from RBC/UKQCD for on-shell K — s- kk matrix elements with m K = 420 MeV and zero 
momentum with and without disconnected contributions [64J. Errors are statistical only. 



7. Summary and Outlook 

Results for the simplest quantities in light quark physics are now in impressive agreement. 
To take one example, many groups now have results for quark masses using dynamical ensem- 
bles, multiple lattice spacings, and improved techniques for computing the renormalization factor, 
and there is agreement at the few percent level, with systematic errors under control. This is an 
important achievement for the lattice. Other quantities that are important for flavor physics show 
similarly impressive agreement, and it has become clear that averages are necessary to maximize 
the impact of these lattice results on constraining new physics. I have reviewed the approach to 
averaging that my collaborators E. Lunghi and R. Van de Water and I have adopted [0], and I have 
presented our updated averages for many quantities involving light quarks. More difficult quanti- 
ties like K — > 7171 matrix elements in the AI = 3/2 channel are now within reach, and preliminary 
results with errors at the 10 — 20% have been presented [86, 80]. The AI = 1/2 channel is more 
difficult but may be attainable in the next few years. 
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